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Abstract 

First protostellar cores are young stellar objects in the earliest evolutionary stage. They are hydrostatic 
objects formed soon after the central portions of star-forming cores become optically thick to dust emission. 
We consider their characteristics in the emitted radiation, and discuss their evolution with increasing mass 
of the cores. Particular attention is paid to detailed radiative and chemical processes in the postshock 
relaxation layer located at the surface of the core, where the majority of radiation is emitted. Most of 
the radiation is originally emitted in the dust continuum in mid-infrared wavelength 10 — 30/.im), which 
reprocessed to far-infrared with ~ 100 — 200jum. Although some fraction (~ 0.1) of the radiation energy is 
emitted in the H2O lines at the accretion shock, most is absorbed and reemitted in the dust continuum in 
the envelope. The H2O lines account for at most ~ 1/100 of the observed luminosity. If a cavity is present in 
the envelope due to outflow or rotation, the dust and H2O line emission in the mid-infrared wavelength from 
the shock can be observed directly, or as a reflection nebula. Among forthcoming observational facillities, 
SPace Infrared telescope for Cosmology and Astrophysics (SPICA) is the most suitable for detecting either 
direct or processed radiation from first-core objects. 
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1. Introduction 

Evolution from the collapse of dense molecular cores 
to the formation of protostars has been discussed the- 
oretically for several decades (e.g., Hayashi 1966; Shu, 
Adams, & Lizano 1987). Advances in radio and infrared 
observations have confirmed many aspects of the evolu- 
tion. The earliest phase of protostellar evolution, i.e., 
the first core, however, still eludes observational confir- 
mation though a theoretical prediction was made almost 
forty years ago (Larson 1969). A first core is a hydro- 
static object that forms after a star-forming core becomes 
optically thick to the dust continuum and the tempera- 
ture evolution becomes adiabatic owing to inefficient ra- 
diative cooling. Their roles in protostellar evolution in- 
clude enabling binary formation and/or driving protostel- 
lar outflow (e.g., Inutsuka 2004). The difficulty in finding 
first-core objects is due to the following reasons. First, 
first-core objects are short-lived. After formation, they 
grow in mass by accretion of the envelope matter and 
the temperature inside increases. When the mass of the 
first core reaches about O.OSM© and simultaneously the 
temperature reaches 2000K, molecular hydrogen begins 
to dissociate. Owing to effective cooling by dissociation, 
the first core begins collapse dynamically. Using a typical 
accretion rate of lO~^M0/yr, the timescale in the first- 
core phase is only 5 x 10*^ yr, which is one or two orders 
of magnitude shorter than the Class phase. This means 
that the expected number of the first-core objects is only 
a few hundredth of Class objects. Besides, the fact that 
first-core objects are deeply embedded in their parental 
molecular cores prevents detection. 

Such difficulties are being overcome with a rapidly in- 



creasing number of observed young stellar objects, as well 
as with dramatical improvement of the sensitivity in the 
infrared/submillimeter wavelengths. Encouraged by these 
facts, the search for protostars in their earliest phase is be- 
ing undertaken vigorously. Among them, Onishi, Mizuno, 
& Fukui (1999) found L1527F (also known as MC27) in 
Taurus by the H^^CO"'" emission, which is a candidate 
for an object in the earliest stage of a protostar owing 
to its very high density (~ lO^cm^"^). Recently, a cen- 
tral submillimeter source (L1521F-IRS) and near-mid in- 
frared bipolar reflection nebulae with a luminosity 0.05Lq 
was detected in this object by Spitzer Space Telescope 
(Bourke et al. 2006). Another more promissing object is 
Cha-MMSl, whose evolutionary stage appears to lie be- 
tween L1527F and a starless core L1544 from the degree 
of N2II+ fractionation and from the detection of a mid- 
infrared source that is several times less luminous than 
L1521F-IRS (Belloch et al. 2006). 

Knowledge of expected observational appearence of 
first-core objects is crucial in finding such objects and 
confirming their identity. Some efforts have already been 
made along this line. Boss & Yorke (1995) presented 
the expected spectral energy distribution of an object 
in the first-core phase, using the result of a three di- 
mensional simulation of star formation. By performing 
a one-dimensional radiative hydrodynamical calculation, 
Masunaga, Miyama & Inutsuka (1998) followed the for- 
mation and evolution of the first cores, and presented the 
evolution of their emission spectrum. The above men- 
tioned hydrodynamical studies, however, failed to follow 
the structure of the postshock flow of the accretion shock 
due to its extreme thinness, although the majority of emis- 
sion is created there (Stabler, Shu, & Taam 1980). In ad- 
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dition to this, those authors included only the dust ther- 
mal emission in their calculation. In this paper, we take 
another approach. With a simpler assumption concern- 
ing the dynamics of the accreting envelope, we correctly 
solve the structure and emission process in the postshock 
flow. Using the injection radiation from the shock, we 
solve radiative transfer in the protostellar envelope. In 
this way, we obtain the observable radiation spectrum 
from the first-core object. 

The organization of this paper is as follows. In Section 
2, we describe the model for calculating the emission from 
the first-core objects. In Section 3, we present the results. 
In Section 4, we summarize this paper, and present a dis- 
cussion and the conclusion. 

2. Model 

In Figure 1, we present a schematic view of our model 
of a first core and its envelope. For the structure and 
evolution of the envelope, we adopt a self-similar solution. 
Given the radius of the first core, which is the location 
of the accretion shock, the shock velocity and preshock 
density are obtained from a self-similar solution. Chemical 
and radiative processes in the postshock relaxation layer 
are analyzed by assuming a steady state shock. Finally, 
using the radiative energy injection obtained as mentioned 
above, we solve the radiative transfer in the envelope. We 
describe specific procedures in the following. 

2.1. Envelope density distribution 

We adopt a similarity solution as an envelope den- 
sity structure. Two types of similarity solutions, (i) the 
Larson-Penston (LP; Larson 1969; Penston 1969) and (ii) 
Shu (Shu 1977) solutions, are studied. In translating 
the similarity solutions to physical quantities, we assume 
isothermal evolution at lOK. In Figure 2, we show the 
density distributions in the envelope both for LP and Shu 
similarity solutions. The structure at three epochs, where 
the mass of the first core Mpc = 0.0125,0.025 and O.O5M0, 
is illustrated for each model. As the maximum mass of the 
first core, we take Mpcmax = 0.05AfQ following Masunaga 
& Inutsuka (2000). As the total mass of the core, i.e. 
the sum of the first-core mass and the envelope mass, 
we take IMq. In both solutions, the density low follows 
p oc. r~^^^ for inner radius and becomes steeper, (x r^^, 
outside. The rarefaction wave, which separates the inner 
and outer parts, propagates outward at the sound speed. 
In the outer region, the Shu solution coincides with the 
singular isothermal sphere, while the density of the LP 
solution is 4.4 times higher. 

Note that even with the same Mpc, the time which 
elapsed since the formation of the first core is differ- 
ent for the two solutions. Since the mass- accretion rate 
is M = 46.9cl/G = 7.5 x lO-^Mg/yr (M = 0.975c3/G = 
1.6 X lO~^M0/yr) for the LP (Shu, respectively) solution, 
which is constant for isothermal similarity solutions, the 
age of the first core is tpc ~ Mpc/M, and its entire life- 
time is tpcmax = Mpc.mnK/ M = 670yr for the LP solu- 
tion, and 3.2 x lO'^yr for the Shu solution. For the same 



first-core mass, the rarefaction wave propagates to about 
a 50-times larger radius in the Shu solution than in the 
LP one. Thus, the LP solution is hardly modified, except 
at the innermost region, while the Shu solution evolves 
remarkably during the first-core stage (see Figure 2). 

We take the radius of the first core, i?sh, as a free pa- 
rameter. For the fiducial case, we take i?sh = 5AU from 
the spherical symmetric model (Masunaga, Miyama, & 
Inutsuka 1998). Wc fix the radius in time because it is 
found to be almost constant from hydrodynamical simu- 
lations (Masunaga ct al. 1998). Taking account of the 
possibility that rotating first cores can be more fiattcned 
and extended (Saigo & Tomisaka 2006), we also study the 
case of Rsh = lOAU in order to mimic this effect. We then 
consider four models that are combinations of either LP 
or Shu solutions, and either R^h =5, or 10 AU. These four 
models are called as LP05, LPIO, SOS, and SIO, hereafter. 
At a given time, the mass inside Rsh is used as the first 
core mass, Afpc- 

2.2. Structure and radiative processes in the postshock 
relaxation layer 

The postshock flow is treated as a steady-state shock 
problem. The preshock density is given by a self-similar 
solution once the shock radius is fixed. The preshock tem- 
perature is given by solving radiative transfer in the en- 
velope (Sec. 2.3). 

Using the shock velocity Vs, Mach number Ai, preshock 
pressure pi and number density ni, the postshock initial 
values can be obtained by the following Rankin-Hugoniot 
relations: 

P2 ^ 27adA^^ - (7ad - 1) 

Pi 7ad + l ' ^ ' 

V2^rH^ (7ad-l) + 2/A^^ 

Vs n2 7ad + 1 ' 

where we use subscript 1 (2) for the preshock (postshock, 
respectively) variables, and 7ad is the ratio of specific heat. 

In the postshock relaxation layer, since the gas kinetic 
temperature T is higher than the radiation temperature 
Trad I the thermal energy is converted to the radiation en- 
ergy. We follow the evolution of a postshock gas parcel 
until T reaches l.lTrad in the following way. The equa- 
tions of conservation of mass and momentum fluxes, 



PV = P2V2=J, 

pv'^ + p = P2vl+p2, 



(3) 
(4) 



can be incorporated in a single equation, 
f 

p+ — = P2V2 +P2 = const. (5) 
P 

Coupling the above relation with the equation of state for 
an ideal gas 

pkBT 



P = 

pma 

we solve the energy equation 



(6) 



No. 



Observational Characteristics of the First Protostellar Cores 



de d 1 1 



(7) 



(8) 



where the energy density per unit mass is 
kBT 

6 = 

(7ad - 1)mtoh ' 

and A is the radiative cooling rate. This determines the 
structure of the postshock flow. 

The radiative cooling is by line emission and dust ther- 
mal emission, 



A = Aline + A„ 



(9) 



In the first term Aiino, we include the fine-structure and 
metastable transitions of C, 0+ and O, and molecular ro- 
tational transitions of H2, CO, OH and H2O. Those are 
treated as in Omukai et al. (2005), but the related coeffi- 
cients for CO and H2O are updated using those in Leiden 
Atomic and Molecular Database (Schoier et al. 2005). 
The cooling rate per unit volume by a transition i is 



(10) 



where n(i„) is the population in the upper level iu, and 
Ai is the spontaneous radiative decay rate. The effect of 
photon-trapping is accounted for by the escape probability 
(Neufeld & Kaufman 1993), which is given by using the 
optical depth at the line centroid Ti 



1 



l + 3r, 



and 



Ti{s) = / ai{s')ds', 



(11) 



(12) 



where s is the distance from the shock. The above intcrgal 
is over flows with \v{s') ~v{s)\ < [Au£i(.s') + AiJ£)(s)]/2, 
where Avd is the thermal Doppler width. The second 
term in equation (9) is the cooling rate by the dust con- 
tinuum 



Agr ^AtT j [S^(Tdust) - By{T^^i)] dv, 



where is the dust opacity and Tdust is the dust tempera- 
ture. The dust opacity is taken from the IPS ("iron-poor" 
silicate) composite aggregates model of Semenov et al. 
(2003), who extended Pollack et al. (1994)'s calculation. 
The dust temperature is calculated by the energy balance 
between radiative loss from the dust and coUisional heat- 
ing by gas (see Omukai et al. 2005 for details). The 
postshock flow is optically thin to the dust continuum. 
Among radiative cooling processes mentioned above, the 
H2O emission and dust thermal emission arc dominant. 

Chemical processes among compounds of H, C and O 
arc included, as in Omukai (2000), where the reaction rate 
coefficients for the C and O chemistry are based on Millar 
et al. (1997). Hydrogen is assumed to initially be fully 
molecular. The initial abundances of other species in the 
gas phase arc set to be j/(CO) ~ 1.10 x 10~^, j/(C) = 1 x 
10"*^, y(0) ^ 3.31 X 10"", and y(H20) = 1.7 x 10"^ follow- 
ing Pollack et al. (1994) 's "Cloud Cores" model in their 
Table ID. If the postshock temperature is higher than 



the evaporation temperature of ice, which is about 150K, 
molecules locked in the ice are released into the gas phase. 
In this case, the amounts of H2O, CO and CO2 in the ice 
mantle ?/icc(H20) = 9.3 x 10-^ yice(CO) = 2.8 x 10-^ and 
yicc(C02) = 5.9 x 10~^ are added to the gas-phase chem- 
istry. 

2.3. Radiative Transfer in the Envelope 

To solve the temperature structure in the envelope and 
the final processed radiation, we need to specify the injec- 
tion radiation from the accretion shock, which is the sum 
of the emission from the postshock relaxation layer and 
radiation from interior of the first core. In the previous 
section, we described the method to calculate the emission 
from the relaxation layer. At the bottom of this layer, the 
radiation and matter are thermally well-coupled. Thus 
the radiation intensity is given by the blackbody radia- 
tion at temperature Trad = Tpost , where Tpost is the tem- 
perature there. To find the post relaxation temperature 
Tpost, we follow Stabler, Shu, & Taam (1981). For the LP 
solution, the envelope is optically thick both to the shock 
emission and to dust re-emission in the envelope. In this 
case, both the preshock temperature Tpre and the temper- 
ature at the bottom of the postshock relaxation layer Tpost 
are equal to the radiation temperature in the postshock 
layer, which is optically thin. Thus, the condition 



T 



post 



T 

' J- pre 



(14) 



needs to be satisfied. For the Shu solution, the envelope is 
optically thick to high energy shock emission, but optically 
thin to dust re-emission. In this case, the condition for 



Tpost becomes 



T 



post 



post ) 



1/4 



(15) 



where Tpre (Tpost) is the radiation energy fiux in the 
preshock flow (at the bottom of the relaxation layer, re- 
spectively) . 

(13) Using the postshock relaxation temperature Tpost ob- 



tained above, the injection radiation to the envelope from 
the accretion shock is 

Iu;m = ly.stiock + Iu,corc ~ -^u.AxIt^ + -Biy(Tpost), (16) 

where £^^_sh is the energy emission rate in the postshock 
layer per unit area, and can be obtained by integrating 
the radiative cooling rate along the postshock flow. 



Ty.sh = / A^ds = / kyvdt 



(17) 



where is the radiative cooling rate per unit frequency. 
The factor 1/2 in the first term on the right hand side of 
equation (16) comes from the fact that the other half of 
the radiation goes into the first core. 

The radiation field in the envelope is determined by 
solving the radiative tranfer under the input radiation 
from the inner boundary. We used the variable Eddington 
factor method (Mihalas & Mihalas 1984). The tempera- 
ture T is determined by the radiative equilibrium 



K. Omukai 



[Vol. 



',dust 



I 



iJh 



(^j^g) Since the flow velocity is approximately given by free-fall 



where Ji/,dust is the monochromatic mean intensity in the 
continuum, and J; is the mean intensity integrated over 
frequency in a line Since the gas and dust are thermally 
coupled in the envelope, we do not distinguish their tem- 
peratures. In solving radiative transfer, we only consider 
the dust continuum as a source of opacity in the enve- 
lope. In all calculated mass ranges in the LP05 model 
and some in the LP 10 model, the dust ice mantle subli- 
mates and II2O vapor exists in the innermost portion of 
the core. However, this part is also optically very thick 
to dust continuum, and the II2 0-line luminosity is signifi- 
cantly attenuated by the dust absorption in the surround- 
ing region. Since we are not interested in such weak line 
emission, we do not include the II2 0-line opacity in our 
model. 



2.. 



Method of Model Construction 



The model is constructed as follows: 

(i) Given the model parameters, namely, the shock ra- 
dius i?sh and the type of similarity solution (LP or Shu). 

(ii) Calculate the shock velocity Vs and the preshock 
density pi for a given i?sh (§2.1). Then, assume the 
preshock temperature Tpre and the radiation temperature 
in the postshock layer Trad- 

(iii) Solve the postshock flow (§2.2) and obtain the in- 
jection radiation field Iv;nx from the shock. 

(iv) Solve radiative transfer in the envelope (§2.3) and 
obtain T-p^^. 

(v) If the assumed and calculated Tprc coincide and the 
shock boundary condition (eq. 14 or 15) is satisfied, the 
model is constructed. If not, go back to (ii) with modi- 
fied Tpre and Trad, and repeat this procedure until conver- 
gence. 

3. Results 

3.1. Emission in the postshock relaxation layer 

In Figure 3, we present the temperature distribution in 
the postshock relaxation layer for three epochs with the 
first core mass being 0.0125, 0.025, and 0.05 Mq (i.e., 
0.25, 0.5, and 1.0 Mpcmax, respectively). The distribu- 
tion is shown as a function of the column density from 
the accretion shock: a higher column density corresponds 
to more downstream. The cooling rate by each cooling 
agent, namely, the dust thermal emission and the II2O 
and Hj^O line emission, is shown individually in Figure 4 
for the same fiows. 

As can be seen in Figure 3, the postshock temperature 
T2 is more dependent on the assumed first-core radius i?sh, 
rather than the type of similarity solution. It is 700 — 900K 
for the cases of i?sh = 5 AU, while T2 ~ 400K for those 
of i?sh = 10 AU. This can be understood as follows. By 
assuming pivl ^ pi in the preshock and P2V2 in the 
postshock, 
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_7ad_P2 _ 7ad fcBT2 

7ad - 1 P2 7ad - 1 pmn ' 



(19) 



1 2 GMfc 
— V — . 

2 ^ Rsh 



(20) 



the postshock temperature can be written as 



7ad 



O.O5M0y I 5AU/ 



(21) 



Since the equation above does not contain the type of 
similarity solution, the values of postshock temperature 
are similar for the two solutions with the same i?pc and 
Mpc- On the other hand, the dust temperature, which is 
close to the radiation temperature, is higher for the LP 
solutions, because of higher luminosity and higher optical 
depth in the envelope. Assuming the radiative diffusion 
equation 

9T ^ 3 KnpL ^22) 

dr Idirac r'^T^ 
in the envelope and approximating the left hand side by 
Ti/i?onv, where i?onv is the outer radius of the envelope, 
we obtain 

(23) 



' cnv 



where Tonv is the envelope optical depth. This shows that 
the higher luminosity and higher envelope optical depth 
result in a higher radiation temperature (~ Ti) at the 
first-core surface. 

Despite the smaller cooling rates (see Fig. 4), the 
shocked gas cools at a lower column density in the Shu 
models (< lO^Ocm-^) than in the LP ones (> lO^Ocm-^). 
This is because of the longer timescale in the Shu cases: 
since the accretion rates of the LP and Shu solutions are 
different by a factor of about fifty, for the same column 
density the time which elapsed after the shock is also dif- 
ferent by the same factor. 

In all models, for low column densities, thermal en- 
ergy in the postshock flow is dissipated mostly in dust 
continuum emission (see Fig. 4). The postshock tem- 
perature is several hundred degree K (Fig. 3), and in- 
creases with the first-core mass, reaching 900K (LP05) 
for Afpcmax = O.O5M0. In this temperature range, the 
most efficient molecular coolant is II2O, and its cooling 
rate is as large as that by dust continuum for small col- 
umn densities (see Fig. 4). However, for Nh ^ lO^^cm"^, 
the H2O lines become optically thick and the cooling rate 
decreases. The optical depth effect is apparent when we 
compare the H2O cooling rate with that by its isotope 
Hj^O, which remains optically thin. 

In Figure 5, we show the contribution to the luminosity 
emitted in the postshock layer by the individual compo- 
nents; dust continuum, H2O lines, H2®0 lines, and OH 
lines. Although the contribution of the H20-line cool- 
ing to the shock dissipation increases with the first-core 
mass, the maximum fraction of radiated energy in the H2O 
lines is at most about 10%, except for the S05 model. 
In S05 model, the H20-line cooling is as important as 
the dust continuum at the end of the first-core phase 



No. ] 

Mpc = O.O5A/0. In this model, since the postshock gas 
temperature is high (TOOK) and it takes a long time for 
the layer to become optically thick, the H2O lines con- 
tribute to cooling as much as the dust emission. 

The contribution of the water isotope H2^0 is very 
small owing to its small abundance. The second impor- 
tant species as a coolant is OH. However, its role is only 
temporary: since OH is intermediate species in chemical 
reactions leading to H2O formation, its cooling rate de- 
creases as OH is being converted to H2O. 

The value of the H2O luminosity is dependent on the 
H2O fraction as well as the gas temperature. In Figure 
6, we show the H2O and O concentrations at the bottom 
of the postshock relaxation layer. In LP05 model, the 
postshock dust temperature is higher than the evapora- 
tion temperature of the ice, even in our minimum mass 
case of O.O125M0. Thus a large amount of water vapor 
[ y(H20) = 1.1 X lO""* ] is already present just below the 
shock. For the lowest mass case Mpc = O.O125M0, the 
final H2O concentration remains almost the same as the 
initial value (Fig. 6 a). In the cases with higher mass, the 
final H2O concentration increases with Mpc by a chemical 
process in the postshock flow owing to higher temperature 
and thus higher H2O production rate (Wagner & Graff 
1987; Kaufman & Neufeld 1996). For Mpc > O.O25M0, 
most of the remaining oxygen, which is not locked to CO, 
is converted to H2O via OH. This can be observed as a 
decrease in the OH luminosity (Fig. 5 a). In the LPIO 
model, the ice mantle evaprates for Mpc > O.O28A'/0. This 
can be seen as the discontinuous rise in the H2O con- 
centration, as well as the H2O luminosity (Fig. 5 b). 
Additional H2O production by chemical reactions slightly 
increases the H2O concentration in the high-Mpc cases. 
The ice mantles do not evaporate in cither the S05 or SIO 
model. In the SOS model, additional H2O is produced by 
high-temperature gas-phase chemistry for Mpc ^ O.O3M0, 
while in the SIO model H2O hardly increases from the ini- 
tial value of 2/(H20) = 1.7 x 10"^ 

In Table 1, we list the important transitions, their en- 
ergy fluxes inline, sh = ^-Eiine.sh and the luminosities 

Liinc = 47ri?2j^Fii„e,sh, (24) 

generated in the postshock relaxation layer. The emit- 
ted luminosities are higher and the wavelengths of the 
strongest emission lines are shorter for LP / small i?sh 
cases than for the Shu / large R^h cases. In general, the 
emitted radiation from the shock is subject to absorption 
and re-emission in the envelope (Sec. 3.2), and thus these 
values does not correspond to the observable line lumi- 
nosities. However, if cavities are present in the envelope 
due to outflow, rotation, etc., the shock emission can come 
out unattenuated. In this case, the emission lines listed 
in Table 1 are observable probably as mid-infrared nebu- 
losity. 

3.2. Processed Radiation through the Envelope 

A first-core object is enshrouded with an optically 
thick envelope, where radiation emitted from the accre- 
tion shock is absorbed and reemitted as dust thermal 
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radiation. The total optical depths of the envelopes for 
each frequency are shown for the LP05 and S05 models in 
Figure 7. The envelope is optically thick at wavelengths 
shorter than 1mm in the LP05 model (LPIO: 600/im, SOS: 
100/im, SIO: 75/xm, respectively). With increasing mass 
of the first core, the inner portion of the envelope is ac- 
creted and the optical depth decreases. This is not a large 
effect compared with the differences between the LP and 
Shu models. 

Figure 8 shows the temperature distribution in the en- 
velope. For the LP models, the dust photosphere locates 
at ~ 30 AU, outside of which the envelope is optically 
thin and reemitted photons escape freely. The temper- 
ature falls approximately as T cx and T oc r~^/^ in 
optically thick and thin regimes, respectively (e.g.. Sec. 
4.2 of Hartmann 1998). Below its vaporization tempera- 
ture (~ ISO K), H2O ice is an important contributor to 
the opacity. Therefore, for ^ 150K, due to evaporation 
of ice mantle of dust, the dust opacity is reduced. This 
causes the flattening of temperature profile in the inner- 
most region observed in the LP models. For the Shu mod- 
els, where the envelope is optically thin, the temperature 
obeys T oc r~^/^ in the entire envelope. 

Figure 9 shows the spectrum of processed dust radia- 
tion, i.e. the radiation escaping from the envelope. The 
emitted dust radiation = AnR^Yi^^^sh in the shock is also 
shown for comparison. However, the peak of the shock 
emission locates at rather short wavelength, around 10- 
100/im, due to high temperature in the postshock flow, 
while that of processed radiation appears at longer wave- 
length, 100/im-lmm, which reflects the lower envelope 
temperature. 

The envelopes of the LP models are optically thick both 
to the shock and the processed radiation, while those of 
the Shu models are optically thick to the emitted radi- 
ation, but optically thin to the re-emission. Owing to 
this difference, the processed spectra of the LP models 
resemble the blackbody radiation (x (x on the 
Rayleigh- Jeans side, while those of the Shu model have 
flatter spectra in low frequencies L^, oc v^l"^ approximately 
(Fig. 9). This flat spectrum is explained by the following 
consideration. The radiation spectrum from optically thin 
medium is 

L, = 4ttJ K,B,dM oc j.3+/3-(3-«)/™ ^ ^1/2^ (25) 

for density and temperature distributions p oc r^" and 
T oc r~™, and an opacity of k^, ozv^. In the above, we 
used n = 3/2, to = 1/3 and (3 — 2. Comparing the LP and 
Shu models with the same Rsh, we notice that, although 
the peak frequency of the shock emission is higher in the 
LP models, that of the processed radiation is lower owing 
to the larger optical depth of the envelope. 

In Figures 11-14, we show the emitted and processed lu- 
minosities of the H2O lines in three different epochs for all 
the models. For the emitted luminosities, i.e. the unat- 
tenuated luminosities before processing in the envelope 
(eq. 24), the corresponding values of energy flux at the 
shock Fiinc,sh are also indicated on the left axis. We also 
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present in Table 2 the luminosities and the line flux at 
the distance of 150pc for the ten strongest H2O emission 
lines at the epoch of Mpc = O.OSMq. Among H2O emis- 
sion lines from the shock in the LP models (see Figures 
11 and 12), those with frequencies ^ 10^'^Hz (30/im) have 
the highest luminosities ~ 10~^ — 1O~^L0. These lines 
are, however, heavily attenuated in the envelope owing to 
the high optical depth in these frequencies. On the other 
hand, the peak in processed line radiation appears around 
lO^^Hz (300^m), where the optical depth of the envelope 
is about a few. With such small values of the processed 
luminosities 10^^ — 10~^'^Lq, those lines are completely 
buried in the dust continuum: iuno/AuD «Cii/,dust, where 
Awd is the Doppler width of the line. In the Shu models, 
line emission in the shock peaks around 10^'^Hz (30//m), 
which is similar to that in the LP model. The processed 
line emission, however, peaks around 3 x lO^^Hz (100/^m), 
higher than in the LP models owing to the smaller optical 
depth in the Shu models. With maximum luminosities of 
10~^ — 10~^Lo, the monochromatic luminosities of these 
lines are higher than that of the dust continuum. The fi- 
nal processed luminosities in the dust continuum and the 
H2O lines are shown in Figure 10. In the LP models, the 
H2O luminosities are only ~ 10^^ of the dust luminosities, 
while in the Shu models, about 1% of the radiation is still 
in H2O lines after reprocessing in the envelope. 

Profiles of brightness temperature are shown for differ- 
ent wavelength bands in Figures (15)-(18). These figures 
show that in low frequencies where the envelopes are op- 
tically thin, the emission from the first cores is directly 
visible. Then, a discontinuity appears in the brightness 
temperature distribution at the radius of the first core. 
On the other hand, at the shorter wavelengths, where the 
envelope is optically thick, only the emission from the dust 
photosphere in the envelope is visible. With increasing ob- 
serving frequency, the brightness temperature at the cen- 
ter decreases and the discontinuity at the first-core surface 
disappears. 

4. Summary and Discussion 

We calculated the radiation spectrum of young stel- 
lar objects in their first-core phase. As envelope mod- 
els, we studied cases of the Larson-Penston (LP) and the 
Shu self-similar solutions. In the LP model, the emit- 
ted spectrum has a peak at ^ 10/im both for the dust 
continuum and H2O lines. The optical depth of the en- 
velope at this frequency is 100— 1000, such that virtually 
all of the shock emission is absorbed in the envelope and 
reemitted as the dust thermal emission at lower frequen- 
cies. The amount of total luminosity remains unchanged, 
which is 0.01 — O.IL©, but the emission peak is shifted 
to lower frequency owing to the lower temperature in the 
envelope. The processed spectrum can be approximated 
by the blackbody radiation with temperature 30K, which 
peaks at 200/im. Our result for the continuum spectrum is 
consistent with a previous study by Masunaga & Inutsuka 
(2000). In the Shu models, the amount of total luminosity 
is 10^^ — 1O~'^L0, which is smaller by almost two orders of 



magnitude than in the LP model owing to the smaller ac- 
cretion rate by a similar factor. The dust temperature in 
the postshock layer is lower (< lOOK), and thus the wave- 
length of the emission peak is longer AOfim) than that 
in the LP case (~ 10/xm). However, owing to the lower 
optical depth in the envelope, the peak of the processed 
spectrum is at a shorter wavelength (100 — 200/xm) than 
that in the LP case. 

Strong H2O emission from accreting protostars has al- 
ready been pointed out by Ncufcld & HoUcnbach (1994), 
who considered the shock on the surface of circumstcl- 
lar disks of more evolved objects. They carried out a 
parameter study for shocks with a preshock density of 
ni = 10^'^ — lO^^cm^^ and a velocity Vs = 5 — lOOkm/s, and 
calculated the emission from them. Although, in their cal- 
culations, protostars tend to be more massive, and shocks 
are thus stronger than ours {ni ~ 10^ — 2 x lO^-'^cm"'^ and 
Us — 1 — 4km/s in our models), our shock calculations 
and results are very similar to their weakest shock cases. 
Provided that shock conditions are similar, it is a difficult 
task to distinguish observationally whether the shock is on 
the disk surface or on the stellar surface without spatially 
resolving the emission region. In the case of accretion onto 
the disk, the amont of dust extinction can be smaller than 
that in our case. Using this fact, we might be able to tell 
the geometry of the shock. We need more studies in the 
future to confirm this possibility. 

So far, we have discussed the case of spherical symme- 
try, where the envelope absorbs the shock emission and 
reemits it at the longer wavelengths. In practice, however, 
outflows are observed almost ubiquitously in very young 
stellar objects, such as Class objects. Theoretically, it is 
known that outflows can be launched as early as the first- 
core phase (Machida et al. 2006). Even without outflow, 
a flattened disk-like density distribution is realized by the 
presence of angular momentum (Terebey, Shu, & Cassen 
1984). Thus, it is highly plausible that the envelope is 
optically thin in the polar direction while it is optically 
thick in the equatorial (disk) one. If so, the shock emis- 
sion escapes unattenuated and can be observed directly 
in the polar direction, or as a reflection nebula in other 
directions. In these cases, the radiation is mostly emitted 
in the mid-infrared range, and about 10% of the emission 
is in II2O lines at similar wavelengths (see Table 1). 

In this paper, we demonstrated that observations in 
mid- to far-infrared and submillimeter wavelengths are 
important for finding the first-core objects. Among 
observational facilities at these wavelengths, we ex- 
amine observational feasibilities by currently operat- 
ing Spitzer Space Telescope (Spitzer) and Stratospheric 
Observatory for Infrared Astronomy (SOFIA), as well 
as forthcoming Herschel Space Observatory, Atacama 
Large Millimeter / submillimeter Array (ALMA) and 
SPace Infrared telescope for Cosmology and Astrophysics 
(SPICA). Spitzer has a spectroscopic line sensitivity of ~ 
10~^^W/m^ in mid-infrared wavelengths (~ 30/im) and a 
photometric sensitivity of O.lmJy at 10/im and ^ lOmJy 
at 100/im. Here we suppose the observation time to be one 
hour. The airborne telescope SOFIA has a spectroscopic 
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line sensitivity of ^ lO^^^W/m^ at mid to far-infrared 
wavelengths (5 — 200//ni) and photometric sensitivity of 
~ ImJy at 10/im and ~ 30mJy at 100/^m. Herschel is a 
space telescope with its aperture 3.5m and oberving wave- 
lengths 60 — 700/im. Its spectroscopic line sensitivity is 
10-1* - 10-"W/m2 in short wavelengths (A < 200^m) 
and 10"^^ — lO^i^W/m^ in longer wavelengths, while its 
photometric sensitivity is 1 — 10m Jy. The millimeter and 
submillimeter telescope ALMA observes at 300^m — 1cm 
(or 30 — 950GHz), with its spectroscopic line sensitiv- 
ity ~ IQ-^^W/m^ and photospheric sensitivity ^ O.lmJy. 
Finally, SPICA is an infrared space telescope with a 3.5m 
aperture, the same as Herschel, but its mirror is cooled to 
4.5K, far cooler than that of Hersehel (80K). The observ- 
ing wavelengths are from 5/im to 1mm, with its spectro- 
scopic line sensitivity being 10^^^ — lO^^W/m^. Its pho- 
tospheric sensitivity depends strongly on the wavelength: 
it is ~ l^^Jy around lO^m and ~lmJy around lOO/im. 

Here, we evaluate the H2O line flux from first-core ob- 
jects. After attenuation in the envelope, the maximum 
value of the H20-line luminosities is iunc ~ lO-^i© with 
its wavelength around 100/im in the Shu models, although 
it is significantly dimmer in the LP models. If there is a 
cavity in the envelope owing to outflow or rotation, and 
the shock radiation leaks through it without reprocessed 
in the envelope, the line luminosities take the maximum 
value at around 30/xm both for the LP and Shu solutions, 
with its value being Lunc ^ \Q~^Lq (LP) or Liinc ^ lO^^L© 
(Shu) , respectively (see Figures 11 - 14). Below, as a 
typical distance to a first-core object, we use Ikpc, far- 
ther than neighboring star-forming regions, because, for 
detection of such rare objects, the observation of a large 
number of young stellar objects is needed. The line flux 
can be given by 

fl„„.3,4xl„-.w(j^)(j|;)" (26) 

= 3.4xlO-W/„^(j^)(jiij)".(27) 

where we use typical values of the luminosities for unat- 
tenuated (IQ-^Lq) and attenuated radiation (lO^^L©) in 
the first and second expressions, respectively. If the shock 
emission from a first core is visible without attenuation in 
the envelope, Spitzer is able to detect it in sources as far 
as Ikpc. If sources are closer, SOFIA also has chance to 
observe them. However, their sensitivities at far-infrared 
wavelengths are not enough for observing the attenuated 
emission. On the other hand, SPICA can detect both 
processed line radiation through the envelope and direct 
radiation from the accretion shock as far as Ikpc. Owing 
to lack of the sensitivity at the mid-infrared wavelengths, 
the other future telescopes cannot observe shock radia- 
tion. As for the attenuated line radiation, which appears 
in the submillimeter range, the sensitivity of Herschel is 
not sufficient. Although sensitivity of ALMA appears to 
be sufficient, absorption by telluric water vapor prohibits 
the observation of H2O lines in the submillimeter range. 
Next we discuss the observability of the dust contin- 
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uum. Using a typical value vLy = IO^^Lq at a typical 
wavelength of processed radiation A = 200/^m (see Figure 
9), the flux of dust continuum is 

which falls below the sensitivities of currently operating 
SOFIA and Spitzer, but is above those of all three future 
telescopes discussed here. In particular, having sensitivity 
both for direct and processed radiation, SPICA is most 
suitable for observing the first-core objects not only by 
the H2O emission, but also by the dust continuum. 

In summary, SPICA can observe first-core objects as 
far as Ikpc by far-infrared dust emission and (if the col- 
lapse is Shu- like) by the H2O emission. In addition, if the 
hot dust emission and H2O high-excitation-line emission 
in the mid-infrared wavelengths originating from the ac- 
cretion shock come out unattenuated through a possible 
cavity, SPICA can detect them as direct or reflection light, 
depending on the orientation between the cavity and the 
observer. 
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K. Saigo, and Y. Yonekura for suggestions, and the anony- 
mous referee for improving the manuscript. This study is 
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Table 1. Most luminous H2O lines emitted in the postshock layer before processing in the envelope at the epoch of Mpc = O.OSMq 




Fig. 1. First core and its environment (schematic). This shoves our model geometry of the first core and the envelope. We calculated 
the processed radiation through the envelope, in addition to the emission from the accretion shock. 
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Table 2. Most luminous H2O lines after processing in the envelope at the epoch of MpQ = 0.05Mq 




radius r (AU) 



Fig. 2. Number density distributions in the envelope. Depicted are both of the LP and Shu solutions for three epochs of the 
first-core mass Mpc = O.O125M0 (dashed), O.O25M0 (dotted) and 0.05Mq (solid). 
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Fig. 3. Temperature distributions in the postshock relaxation layer. The panels are for (a) LP05, (b) LPIO, (c) S05, (d) SIO models. 
For each model, those for three epochs of M-pc = 0.0125Mq (bottom), O.O25M0 (middle) and O.OSMq (top) are shown. The solid 
curves indicate the gas temperature, while the dashed ones are for the dust temperature. 
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Fig. 4. Cooling rates in the postshock relaxation layer. The panels are for (a) LP05, (b) LPIO, (c) 305, (d) SIO models. For each 
model, three epochs of A/pc = 0.0125A?f© (bottom), 0.025Mq (middle) and 0.05Mq (top) are shown. The solid, dashed, and dotted 
curves indicate the cooling rates by dust, H2O, and Hj^O, respectively. 
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Fig. 5. Emitted luminosities in the postshock layer in the dust continuum (solid), H2O lines (long-dashed), Hj^O lines (short- 
-dashed), and OH lines (dotted) as a function of the first-core mass. The panels are for (a) LP05, (b) LPIO, (c) 305, (d) SIO 
models. 
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Fig. 6. H2O (solid) and O (dashed) concentrations at the bottom of the postshock layer as a function of the first core mass. The 
panels are for (a) LP05 and 10, (b) 305 and 10 models. 
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Fig. 7. Optical depth of the envelope for models LP05 and 305 as a function of the wavelength. The dashed and solid lines show 
the cases of A/po = 0.0125Mq and 0.05Mq, respectively. The dotted horizontal line indicates where the optical depth is unity. 
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Fig. 8. Temperature distributions in the envelope. The dashed and sohd hncs show the distribution at epochs of MpQ = 0.0125Mq 
and O.OSMq, respectively. Lines for T oc r~^^^ and oc r~^, which correspond to optically thick and thin cases, respectively, are also 
shown. 
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Fig. 9. Spectral energy distributions in the dust continuum of the first-core objects. Both the emitted (thin) and processed (thick) 
spectra are shown. The panels (a)-(d) show models LP05, LPIO, SOS, and SIO, respectively. In each panel, the distributions are 
shown for three difl'erent epochs of Afpc = O.O125M0 (dashed), O.O25M0 (dotted), and 0.05Mq (solid). 
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Fig. 10. Evolution of luminosities in the dust continuum (solid) and H2O lines (dashed) at the outer boundary of the envelope, 
shown as a function of the first-core mass. The panels are for (a) LP05 (thick) and LPIO (thin) models, (b) 305 (thick) and SIO 
(thin) models. In panel (a), the luminosities in H2O lines Lh^O multiplied by 10^ for clarity. 



K. Omukai 



[Vol. 



log u (Hz) 
13 12 




X(,um) 



log V (Hz) 
13 12 




A(jU,m) 



log u (Hz) 
13 12 




X(/j,m) 



Fig. 11. Emitted (upper halves) and processed (lower halves) luminosities of H2O lines for the LP05 model. Three epochs of (a) 
O.OI25M0 (b) O.O25M0 and (c) 0.05Mq are shown. In the upper halves, the values of the energy flux at the shock are also indicated. 
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Fig. 12. Same as Figure 11, but for the LPIO model. 
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Fig. 13. Same as Figure 11, but for the SOS model. 
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Fig. 14. Same as Figure 11, but for the SIO model. 
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Fig. 15. Distribution of the brightness temperature at lOOGHz as a function of the impact parameter of the Hne of sight. Panels 
(a)-(d) show LP05, LPIO, SOS, and SIO models, respectively. In each panel, the distributions are shown in three different epochs 
Mpc = O.O125M0 (dashed), O.O25M0 (dotted), and O.OSMq (solid). 
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Fig. 16. The same as Figure 15 but for the observing frequency at 300GHz. 
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Fig. 17. The same as Figure 15 but for the observing frequency at 860GHz. 
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Fig. 18. The same as Figure 15 but for the observing wavelength at 100/im (frequency at 3THz). 



